Plasma instabilities as a result of charge exchange in the 
downstream region of supernova remnant shocks 
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ABSTRACT 

Ha emission from supernova remnants (SNRs) implies the existence of neutral 
hydrogen in the circumstellar medium. Some of the neutral particles penetrating 
the shock are ionized by the charge exchange process and make a cold ion beam in 
the shock downstream region. We perform linear analyses of collisionless plasma 
instabilities between the cold beam and the hot downstream plasma. We find 
that, under typical SNR conditions, either the resonant instability or the Weibel 
instability is the most unstable. This mechanism may amplify the magnetic field 
to more than 100/iG and changes the shock structure. As a result, the radio 
spectrum and the large magnetic field can be explained, apart from the widely 
discussed Bell's mechanism. 
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1. Introduction 



Supernova remnants (SNRs) are observed in various bands from radio to gamma-ray. 
X-ray observations provi de the evidence that electrons are accelerated to highly relativistic 
energies in SNR shocks (iKoyama et al.lll995l ). Moreover, it is believed that ions are also 
accelerated to the "knee" (~ 10^^ eV) of the Galactic cosmic ray spectrum as the origin 
of galactic cosmic rays. X-ray observations of SNRs also suggest the ma gnetic field am- 
plific a tion up to 10^~^/xG i n the vicinity of the collisionless shock fr ont (jVink fc Laming 
20031 : iBerezhko et al.l l2003l : iBamba et al.l l2005l : lUchiyama et al.l 120071 ) . Theoretically, sev- 
eral mechanisms for magnetic field amplification have been proposed, one of which is the 
cosmic ray driven instability, whereby cosmic rays diffuse to the upstream region where the 
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resonant and the nonresonant instabilities become unstable and cause fluctuations in the 



magnetic field flLucek fc Belli I2OOOI : Belli 12004 : iReville et al.ll200d 120071 ). Rec ently, the non- 



resonant instability has been i nvest i gated with Particl e In Cell simulations (iNiemiec et al. 



2008 



2005 



Reyille et al. 



Riauelme fc SpitkovskyI boosi: lohira et al.l [2009[ ) and MHD simulations dSeU 



2004 



20081 : Izirakashvih fc PtuskinI l2008h . but lohira et al.1 J2009h claim that 



the periodic boundary simulations cannot describe nonlinear features. Thus, the ampli- 
tude and the wave spectrum of the magnetic fields excited by the nonresonant instability in 
the downstream region are open issues. The other r nagnetic field amplification mechanism 



is fiu id type instabilities in the downstream region (jGiacalone fc Jokipiil l2007l : llnoue et al. 



2OO9I ). These instabilities generate waves with longer wavelength than plasma kinetic insta- 
bilities. 

In this letter, we propose yet another mechanism of the magnetic field amplification, 
where neutral particles in the partially-ionized interstellar matter (ISM) play a significant 
role. The existence of neutral particles has been identified in many young SNRs from the 
detection of Ha emissions. Since SNR shocks are collisionless, neutral particles can pene- 
trate the shock front freely without deceleration. As a result, a significant relative velocity 
exists between the plasma particles and neutral particles in the shock downstream region, 
which is given by frei=3vsh/4 {vgh is the shock velocity). Through the charge exchange pro- 
cess in the downstream region, shocked hot protons and cold neutral hy drogen atoms are 
converted to hot neutral hydro gen atoms and cold protons, respectively (IHeng fc McCray 
20071 : Ivan Adelsberg et al.ll2008l ). It is noted th at the Ha emissions from SNR consist of two 



components with narro w and wide line widths (j Chevalier fc Raymond 1 19781 : IChevalier et al. 



1980l : ISmith et al.lll99ll ) . The narrow- width component originates from cold hydrogen atoms 
(before charge exhange) excited by heated proton and electron collisions, while the latter 
wide- width component originates from hot hydrogen atoms produced by the charge exchange 
process. 

The neu tral fraction of the ISM aroun d supernovae is often found to be order of unity. 
For example, IChavamian et al.l (I2OOOI . l2002l ) have shown that the neutral fractions are about 
0.8 and 0.1 in Tycho and SN1006 environments. Therefore, the charge-exchanged cold pro- 
tons can have energy density comparable to that of the hot shocked protons, so that their 
effects on the downstream physical condition should be quite important. Such a condition is 
similar to interacting region between the solar wind and a comet or interst ellar neutral par- 
ticles where several collisionless plasma instabilities have been considered (IWu fc Davidson 
1972h . 



We first provide the background plasma condition (§2). We then perform linear analyses 
of the collisionless plasma instabilities for a parallel shock (§3) and discuss impacts on the 
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magnetic field amplification, the shock structure and the particle acceleration (§4). Studies 
of the perpendicular or oblique shocks will be presented in another paper. 



2. Background plasma condition in the shock downstream region 

Although collisionless shock structures have not been understood completely, the dissi- 
pation scale is thought to be about the ion gyroradius, j ~ 10^~^°cni , which is much shorter 
than the ionization scale, kon ~ 10 cm f ivan Adelsberg et al.l[2008l l Therefore the charge 



exchange process occurs after the upstream plasma density is compressed according to the 
Rankin- Hugoniot relatio n. Interestingly, th is ionization scale is comparable to the spatial 



scale of X-ray filaments (IBamba et al.ll2005l ). In this letter, we treat only the parallel shock, 
so that the magnetic field of ISM -Bism is not compressed. For the particle components in the 
downstream region, we consider neutral hydrogen atoms (density nn), shocked hot protons 
(^p,hot) and electrons (rze), and charge-exchanged cold protons (npcoid)- We choose typical 
values for young SNRs as, 

np,hot = ne = Icm"^ BisM = 3yuG, (1) 

and set Vsh = 0.01c. The density of cold protons rip^coid becomes higher as the neutral 
hydrogens flow to the downstream, and is a function of the elapsed time after passing over 
the shock front, 

Dflp cold 

= nnnp^hotCTCEVreh (2) 

where (Tq.e. is the cross section of charge exchange between the hydrogen atom and shocked 
proton. We set ac.E. = 10~^^cm^ as a typical value, so that the charge exchange time 
(^p.hotCTc.E.'i^rei)""'^ = 4.4 X lO^sec. From equation (2), we have 

'^--'^ = ^™ ( 4.4xmec ) ' 

where we set nn to be constant at nno, the atomic hydrogen density in the ISM. This 
approximation holds as long as rzp^coid -C nno- To satisfy the current neutrality, in the 
downstream electron plasma rest frame, drift velocities of two component of protons are 

S'T-pjCold 

^^d,hot — —77 ; r^^sh 

4(np,hot + '^p.coid) 

3^P,h0t /,N 

^d,cold = 77 ■ rVsh- (4) 

4(np,hot + np,cold) 

The temperatures are 
3 

^P,hot = Y^^p^sh = 18.75keV, Tp^coid = leV, and Tg = 0.03Tp,hot, (5) 
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where we assume the downstrea m electron and proton temperature ratio Tp/Tp^hnt = 0-03 as 



the ty pical value of young SNRs (IGhavamian et all2007l : lvan Adelsberg et all2008l : IOhira fc Takahara 
20081 ). 



3. Linear analysis of collsionless plasma 



In this section, we perform linear analysis by treating np com as a free parameter. We 
consider only the cases where the wave vector k is purely parallel or perpendicular to the 
magnetic field. To obtain the maximum growth rate of the instabilities as a function of 
^p,coid/''^p,hot5 we solve numerically the stand ard collsionless p lasma dispersion relation with 
the drift velocity and the thermal velocity (jlchimaru 1 119731 ) because the thermal effect is 
important in the downstream region. 



3.1. k II B 



For parallel propagation, there are three unstable modes. One is the ion acoustic insta- 
bility which is an electrostatic mode. The dispersion relation is given by 



= 1 



1 



uj — kvd. 



-dz. 



(6) 



kv 



th„ 



where the subscript s represents particle species, here electron, hot proton, and cold proton, 
and Ljp^s,Vd,s and Vth,s are the corresponding plasma frequency, drift velocity and thermal 
velocity, respectively, and Z{^s) is the plasma dispersion function. 

The others are the resonant and the nonresona nt instabilies which are electromagnetic 
modes and can amplify the fluctuating magnetic field (IWinske &: Leroyil 19841 ). The dispersion 
relation for right handed electromagnetic waves is given by 







s 

UJ - kVd,s + ^c,s 



UJ 



kv 



th,s 



kv 



th,s 



(7) 
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where Qc,s = qsB/msC is the corresponding plasma cyclotron frequency {q^ — —q^ — e). The 
resonant instability satisfies the condition Re[u}\ — /cwd.cow + ^c,p = 0, so that the propagating 
direction is the same as the drift velocity of cold protons. On the other hand, the nonresonant 
instability does not satisfy the resonance condition and the propagating direction is opposite 
to that of the resonant mode. 



3.2. k±B 



For perpendicular propagation, the Wcibcl instability with the uniform magnetic field 
is excited. In this condition, the unstable mode is the ordinary mode, that is, the fiuctuating 
electric field is parallel to the magnetic field, so that the dispersion relation is given by 







^-(?y-i:(^) 



1 + 2 



n=l 



— ri^^lg 



kv 



th,s 



c,s 



(8) 



where /„ is the modified Bessel function of the first kind with order n. 



3.3. results 

Substituting equations (1), (4) and (5) to the dispersion relations (6), (7) and (8), we 
obtain the growth rates of these instabilities. Figure 1 shows the maximum growth rates of 
these instabilities as a function of np,coid/''^p,hotj where the growth rates are normalized by 
f2c,p- Note that even if 7^p,coid /''^p.hot is large, the growth rate of the nonresonant instability is 
smaller than that of the resonant instability because the thermal velocity of the downstream 
plasma is comparable to the relative velocity 3t;sh/4. The Weibel instability is stabihzed 
when its growth rate becomes smaller than the proton cyclotron frequency because the 
plasma motion is controlled by the uniform magnetic field. In our plasma condition, the 
Weibel instabihty is stabihzed when rip^coid/^'p.hot is smaller than about 3 x 10~^. This critical 
density depends on thermal velocities, the uniform magnetic field and the drift velocities. 

To find which instability becomes nonlinear first after the shock passage, we compare 
the growth rate of the instabihties with the growth rate of the cold proton density t~^. From 
equation (3), we plot l/(iJlc,p) as a function of 7T'p,coid/^p,hot in figure 1. The black solid and 
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dashed lines show the case of who = lcm~^ and 20cm~^, respectively. The black solid line 
and the growth rate of the resonant instability cross at ^p,coid/'^p,hot — 5 x 10~^, so that 
first of all the resonant instability becomes nonlinear after passing over the shock front at 
around 50 times the proton cyclotron time. If the number densit y of hydrogen atoms in 
ISM nnn is ~ 20cm~^ as observed in high density clumps in ISM (IHeiles fc Trolandl 12003 



Inoue fc Inutsukal 120081 ) ■ it is not the resonant instability but the Weibel instability that first 
reaches the nonlinear phase. Even when nno is as high as 20cm~^, the co llision frequency 
with neutral particles is about 10~^sec^^ ~ lO^^^cp (IReville et al.l 120071 ) which is much 
smaller than the growth rates, so that the coUisional damping by the neutral particles is not 
important. Depending on the ISM density, either of the resonant or the Weibel instability 
works to amplify the magnetic field. If this is the case, the high magnetic field intensity found 
in some SNRs is explicable. Although figure 1 shows that the ion acoustic instability is the 
most unstable mode when np^coid/''i'p,hot > 10~^, in fact it is not because cold protons have 
already been slightly heated by the waves excited by the resonant or the Weibel instabillity 
before rip,coid/'"'p,hot attains such a value. For example, when np,coid/''^p,hot becomes 0.1, if the 
temperature of cold protons becomes larger than about 50 times the initial temperature, the 
growth rate of the Weibel instability becomes larger than that of the ion acoustic instability. 
Even during such a large temperature increase of the cold protons, the kinetic energy of 
neutral atoms is almost conserved, so that its further conversion to the magnetic field energy 
remains possible. 



4. Summary and Discussion 

We performed linear analyses of coUisionless plasma instabilities caused by cold proton 
beam generated by charge exchange process. We show that the resonant instability or the 
Weibel instability becomes nonlinear after passing over the shock front and the cold beam 
carries a large kinetic energy. For the estimation of the magnetic field amplified by the 
instabilities, we can us e a simple order-of-magnitude estimation based on the simulations by 



Kato fc Takabd fl2008f ). 



6B ^ 100 ( ^) ( ( /.G, (9) 



0.01/ V0.4cm-3/ VO.Ol 

where is energy conversion parameter from the kinetic energy to the magnetic energy. 
Equation (9) is close to the suggestion by the X-ray observations. To make further discussion, 
however, we need more realistic simulations in which charge exchange process is properly 
taken into account. 



Next, we discuss the modification of shock structure by the existence of neutral atoms. 
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Neutral atoms of ISM can penetrate the shock front without the dissipation and the deceler- 
ation. Although we do not know how the mixing occurs after ionization, we speculate that 
eventually cold ions and hot ions have the same drift velocity which is one of the center of 
mass and that the relaxation scale is the ionizing scale. In this case, the far downstream 
velocity in the shock rest frame fdown is 

1 



^^down = , , „ Vsh, (10) 

where x is the ionization fraction of ISM. Hence, the spectral index of radio synchrotron 
from accelerated electrons a is 



a 



1 
2^' 



(11) 



where the acceleration n iechanism is the diffusive shock acceleration so that a depends on the 
veloc ity ratio Vsh/ Vdown ( iKrymsky 1119771 : lAxford et al.lll977l : lBelllll978l : iBlandford &: Ostriker 
19781 ). If X = 5/6, a = 0.6, it can e xplain the observed radio spectrum slightly steeper than 
the simplest prediction of a = 0.5 ([Reynolds fc Ellisonlll992l ). 



n 



Finally we mention the case where an SNR is in a hot wind bubble with a lower density 
cm~^ and a higher temperature T ~ 10^~^ eV. The gro wth rate of the resonant 
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instability depends on np^coid/''^p,hot and flc,p (jWinske &: Leroylll984j ). The growth rate of the 



1/2 

Weibel instability depends on coid oc n^^^^^ and decreases rapidly when the growth rate 
is comparable to Qc,p- Hence, for lower upstream density, the growth rate of the resonant 
instability does not change from that of figure 1, but the growth rate of the Weibel instability 
becomes smaller than that of figure 1. A higher temperature does not affect much the growth 
rates because the proton thermal velocity is much smaller than the relative velocity between 
cold and hot protons. Thus, the mechanism proposed in this letter works as long as the 
neutral fraction of ISM is as high as ~ 0.1. 
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Fig. 1. — The maximum growth rate of instabihties. The black hne and black dash line show 
l/tQc,p with nno = lcm~^ and nno = 20cm~^, respectively. 



